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We show that the observation of the Drell-Yan production of low-mass lepton-pairs 
(M < 20 GeV) at high rapidities (Y •> 3) at the LHC can make a direct measurement 
of parton distribution functions (PDFs) in the low x region, x < 10~ 4 . We describe 
a procedure that greatly reduces the sensitivity of the predictions to the choice of the 
factorization scale and, in particular, show how, by imposing a cutoff on the transverse 
momentum of the lepton-pair, the data are able to probe PDFs in the important low scale, 
low x domain. We include the effects of the Sudakov suppression factor. 



1 Introduction 

The very high energy of the LHC allows us to probe the parton distribution functions (PDFs) 
of the proton at extremely small x; a region not accessible at previous accelerators. One such 
process is the Drell-Yan production of low-mass, lepton-pairs, at high rapidity pQ , another is C- 
parity-even quarkonia production [2]. Here we study the former process in detail. For example, 
assume, for the moment, Drell-Yan production based on the LO subprocess qq — > 7* — > fi + yT . 
Then the x values of the q, q are 
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for a pair of mass M = 6 GeV and rapidity Y = 4, and collider energy y^i = 8 TeV. Drell- 
Yan events with these kinematics are in reach of experiments at the LHC, see, for example, 
Ref. pj. However, the predictions of the cross section depend sensitively on the choice of 
the factorization and renormalization scales. The dependence on the factorization scale, fip, 
is shown in the left half of Fig. [TJ see also pp. The strong sensitivity is mainly due to the 
dominance of the double logarithmic contribution, (a s ln(l/a;) lnfi 2 F ) n which leads to strong 
scaling violations in the low x domain. Moreover, at a relatively low factorization scale, the LO 
qq — > 7* subprocess is overshadowed by the NLO subprocess gq — > q^*, due to the dominance 
of the gluon PDF at low x. However, it is this very dominance which enables an optimal choice 
of scale to be found pQ which greatly suppresses the scale dependence of the predictions. In 
other words, by choosing an appropriate value of /ip = /i we resum the most dangerous double 
logarithms into the incoming PDFs. In Section [2] we describe the procedure which allows the 
calculation of the value of this optimal scale /io, which should be used in the LO contribution, 
based on the known NLO result. 

In order to probe the parton distributions at lower scales we introduce in Section [3] a cutoff, 
ko, on the transverse momentum, k t , of lepton pair, and study the behaviour of the cross 
section of Drell-Yan events with k t < k Q . Here one faces another problem - large higher-a s - 
order corrections enhanced by logarithms of the ratio ko/M. These double logarithms, of form 
(a s ln 2 (/c /M)) n , can be resummed into the Sudakov T-factor. Thus our aproach accounts 
for (and resums) both possible types of double logarithms which are responsible for the large 
higher-c^-order corrections. As long as the incoming PDFs are known, the remaining NLO 
corrections only weakly affect the theoretical predictions. The cross sections expected from 
these kinematics are evaluated in Section HI Section [5] contains our conclusions. 



2 Procedure to reduce the scale uncertainity 

Basically, the idea in Ref. [I] is to write 

a NLO (/i F ) = PDF( yUF )®(C LO + a s C , r ^ 1 °)®PDF( /UF ), (2) 

and to note that double counting has to be avoided, since part of the NLO contribution is 
already included, to leading log accuracy, in the first-order in a s term in LO DGLAP. The 
effect of varying the scale from m to yUp in both PDFs of the LO-generated contribution can 
be expressed as 

a h °(n F ) = PDF(m) ® (c LO + ^ln O^j (P Mt C LO + C LO P right )) ® PDF(m), (3) 
where the splitting functions 

-fright -Pqq -Pqg > -Pleft Pi qq -Pqg (4) 
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act on the right and left PDFs respectively. We may equally well have incoming g's in P r j g ht and 
incoming g's in P\ c ft- When this LO-generated term is subtracted, the remaining NLO coefficient 
function C^f(/iF), becomes dependent on \ip. As a result, changing \ip redistributes the order 
a s correction between the LO part (PDF <g> C LO <g> PDF) and the NLO part (PDF <g> a s C^? ® 
PDF). 
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Figure 1: (a) Sensitivity of M = 6 GeV Drell-Yan production at the LHC energy of 7 TeV, 

as a function of rapidity Y, to the choice of factorization scale: /xp = M/2, M, 2M, at LO, 
NLO, NNLO, calculated using the programme Vrap [4J, (b) The bold lines correspond to the choice 
fip = fiQ = 1AM which minimizes C^°, and show the stability with respect to the variations 
/i = M/2, M, 2M in the scale of the PDFs convoluted with - the // dependence is indicated 

by the symbolic equation at the top of the diagram. The dashed lines show that the stability 
disappears for other choices of /i The small crosses are the NNLO result. The figure is taken from 

m 



The trick is to choose an appropriate scale, fip = so as to minimize the remaining NLO 
contribution C^°(/ip). To be more precise, we choose a value \ip — fi such that as much as 
possible of the 'real' NLO contribution (which has a 'ladder-like' form and which is strongly 
enhanced by the large value of ln(l/x)) is included in the LO part where all the logarithmically 
enhanced, a s ln(l/x), terms are naturally collected by the incoming parton distributions. 

To implement this for Drell-Yan production at low x, we first note that the majority 

of incoming quarks/ ant iquarks are produced via the low- a; gluon-to-quark splitting, P qg . That 



3 



is, the most important NLO subprocess is gq — y q*y*. The corresponding cross section reads 



dcr(gg — > 37*) 1 a 2 a s 


s t 2M 2 u 


dM 2 dt 9M 2 s 2 


t s st 



a a s z 1 
9M 4 t 



((l- z f + z 2 ) + z< 



M 4 



M 2 



where u = M 2 — s — t ancQz = M 2 /s. The kinematical upper limit of \t\ is 



M 



{l-z)M 2 /z, 



On the other hand, the a s term in the LO-generated DGLAP contribution is 

da(gq -> <n*) _ a 2 a s 2 1 ^ 2 + ^ _ 



(5) 
(6) 

(7) 

(8) 



dM 2 dt 9M 4 t 

We integrate the exact result, ([6]), from |t | up to \ti\, and the approximate expression, Q, 
from 1 1 1 to /Mq, and choose /xo so that the two integrals equal each other; that is 

l + 3z 



z \z 



[z 2 + (1 - z) 2 ] In 



l*o| 



((1 



, M 2 1 
in t— t + In — 



o 



+ 
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where the infrared divergency, as to — >■ 0, cancels. After integration of (|9]) over the incoming 
gluon, quark flux, dz/z, we find that the optimum scale is given by 

' ' no = I ABM. (10) 
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that is 



Indeed, we can check that the explicit calculation of the pure logarithmic LO evolution in- 
tegrated up to fio = 1A5M convoluted with C LO reproduces the NLO contribution of this 
subprocess. Since in the \ow-x region gq — y qj* is the dominant NLO subprocess we anticipate 
that the choice of \ip = /i = 1.45M will minimize the NLO contribution, C^fi^fip). The 
insensitivity of the prediction of the NLO cross section, da/(dMdY), to the choice of scale is 
shown in the right half of Fig. [TJ We see that for fj,Q = 1AM the NLO predictions, using 
the different scales \ip = M/2, M, 2M in the PDFs of the NLO contribution, are virtually 
indistinguishable from each other. 

The predictions use the MSTW2008 NLO set of PDFs j5]. The corresponding PDF uncer- 
tainty for Y > 3 (that is < 4 x 10 -5 ) is such that the error corridor embraces the predictions 
of the other recent sets of PDFs, see the plot in Fig. 5 of [TJ . 

Recall that at lower scales the uncertainties of the present low x PDFs strongly increase. 
Pure DGLAP extrapolations in this small x domain are unreliable due to the absence of ab- 
sorptive, ln(l/x), etc.,.. modifications. Rather, LHC data for this process will provide direct 
measurements of parton distribution functions in this, so far unexplored, low x domain, albeit 
at the moderately high scale \ip — fio — 1.45M. For example, the observations of the produc- 
tion of Drell-Yan pairs of mass M = 6 GeV probe PDFs at a scale Q 2 = /Xq — 75 GeV 2 . 
Clearly, it is desirable to select Drell-Yan data which can measure small x PDFs at lower Q 2 
scales. Can this be done? 



1 Strictly speaking, z is the ratio of the light-cone momentum fraction carried by the 'daughter' quark to that 
carried by the 'parent' gluon, z = x^jx^. 
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3 To probe small x PDFs at low scales 



In this Section we investigate the possibility of probing PDFs at smaller Q 2 scales by imposing 
a kinematical cut on the Drell-Yan events. Naively we may expect that the production of Drell- 
Yan pair with a small transverse momentum of the pair k t < k should be described by 

the PDF measured at scale \i ~ k , since the corresponding distribution includes the partons 
with transverse momentum p t < [i. 

More precisely the leading order k t distribution is described by the formula [6] 
da d 



dM 2 dYdk 2 dk 



[PDF(fc t ) C LO (M) T 2 (k u M) PDF(fc t )] , (11) 



where factor T q reflects the probability not to emit additional gluons with p t > kt during the 
incoming quark evolution up to scale M. This probability is given by the resummation of the 
virtual (oc 5(1 — z)) loop contributions in the DGLAP equation [TJ, [8] 

r fl (*?y) S exp(- f 2 ^^ 1 f ld ^zP qq (z)e(l-z-A)Y (12) 



where the infrared cutoff is 



A = — l —. (13) 
fi + kt 

Note also that selecting the events with a cut on the transverse momentum k t of lepton pair, that 
is, where the heavy photon has k t < ko, we suppress the contribution of the NLO gq — > (77* and 
qq — > (77* subprocesses and therefore increase the relative importance of the LO contribution. 
As a consequence, this will decrease the value we will find for the 'optimal' factorization scale, 
/i . 

As a first step we neglect the role of Sudakov T-factors and calculate the value of fio as a 
function of the cutoff k t > ko, assuming that T = 1. This can be done analytically. Then, in 
the following subsection, we will include the Sudakov effect. 



3.1 Imposing k t (n + n ) < ko reduces the optimal scale //o 

To explore this possibility, we note that, in the framework of the collinear approximation, the 
transverse momentum kt is generated by NLO subprocesses, like qq — > g^* or, most important 
for the low x domain, the subprocess qg — > qY- Repeating the procedure described in Section 
[2j we first calculate the part of the cross section of this last subprocess which satisfies the 
cut k t < k . Then we compare it with the cross section generated by the LO matrix element 
convoluted with the last step of the LO DGLAP evolution up to the factorization scale [ip. We 
choose the value /if = yUo(^o) that, in terms of the LO approach, reproduces the NLO result 
with the k t < k cut imposed. The result of taking the factorization scale \if = ^o(k ) is that, 
as before, we will again absorb into the LO term the major part of the NLO contribution; that 
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is, the part that is enhanced by the large values of ln(l/x). Clearly this 'optimal' scale /xq now 
depends on the transverse momentum cutoff ko, and will decrease if we impose a tighter cut. 

To determine the dependence of /io on the cutoff ko, we must first rewrite the subprocess 
cross section ^ in terms of k\ , rather than t. The relevant relations are 

f 



-t 



s-M 2 ± v /(s - M 2 ) 2 - Askt 



giving 



dt 



Ukj 



dkf. (14) 



Allowing for both the negative and positive roots, (|5| becomes 
d cr(gq — > q-y*) la 2 a s z z 



dM 2 dk 2 



[l + 2z 2 -2z)- 



zM 2 



+ 1 + Sz 



(15) 



9 M 6 ^/(l - zf - Azk 2 /M 2 

We integrate over the k 2 range, equivalent to the t interval t to t\ as before, and then demand 
that the result is equal to the integral over the LO-generated cross section (pO) from |t | to /i . 
We find that the result, equivalent to is 

1 + 3z 



[z 2 + (1 - z) 2 ] In 



l*o| 



((l-,) 2 + , 2 )(ln^+lnl- 



l-z) 



X{z,k Q ] 



(16) 

where again the infrared divergency, as to ^ 0, cancels; and where we now have a new term 

Az k 2 \ 1 + 32 



I 0(z max z^j 



2(1 + 2z 2 - 2z) tanh -1 



(1 



1-2 

M 2 



+ 



/r 2 

; M 2 



(17) 

The optimal scale /i , determined from (16), now depends on the choice of the k t - cutoff, k , due 



to the presence of I(z, ko). What is the origin of the theta function in this new term? Note 
that the fc t -cutoff only has an effect when ko is smaller than the maximum kt- Thus, there is a 
maximum value of z for which the cutoff is to be applied, given by 



1 + 2 



Z-2 

M 2 



-2, 



l- 2 



M 2 \M 2 



k 2 



+ 1 



(18) 



In order to determine the optimal scale /io(&o), we proceed as before, and integrate (16) 
over the parton-parton flux dz/z. Since the gq — > q^* subprocess is described by the same 
diagram as that corresponding to the last step of the LO DGLAP evolution, the scale /i 
that is calculated in this way provides an exact 'resummation' of the most dangerous^] double 
logarithmic terms of the form C n (a s Inq 2 ln(l/x))™ inside the incoming parton distributions. 
After the integration over the flujiQ we find /io(^o) is now given by 

In ( ' l) . (19) 



M 4 V2 



2 That is terms enhanced by the large value of ln(l/x). 

3 In order to account better for the x-dependence of the incoming parton densities we may include in the flux 
an additional factor z s . However, integrating over such a flux factor with 5 < 0.2 gives a negligible change in 
the resulting value of Ho/M for the region of interest fco ^ M/2. 
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in the place of (10). The additional term / can be evaluated analytically. It is given by 



(66w 3 + 72w 2 + 18w + 12) In [(w + l)/w] - 66w 2 - 39w - 4 
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(20) 



(21) 



The dependence of the optimal scale, /i , on the k t cutoff, k t < /cq, is shown by the continuous 
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Figure 2: The dependence of the optimal scale fx on the choice to cutoff k t < k on the transverse 
momentum, k t , of the Drell-Yan pair of mass M. The curves for ko/M > 0.1 show how 

including the Sudakov factor, with and without running a s , changes the dependence of the optimal 
scale fi on the cutoff k . 



curve in Fig. [2j In particular, for a cutoff ko = M/A we find that the optimal scale is 
fio = 0.39M, while for ko = M/2 we have fio = 0.70M. These values are a considerable 
improvement on the result /zo = 1.45M in the absence of the cut, but are somewhat larger 
than the naive estimate of fio = ko since, (i) the scale /j,® is defined in terms of the virtuality 



t 



z) > 



kf, and, 



(n), in the gq — > qY cross section, ([5]), 



\u\ 



but where the transverse momentum k 2 



we have a 'backward 
tu/s is still 



scattering' contribution with \t\ > 
small. 

So it looks as if we have achieved our aim. By selecting data satisfying a cutoff k t (fi + fi~) < 
ko, we do indeed significantly reduce the optimal scale, provided the cutoff satisfies k < M. 
But, first, we must check whether allowing for Sudakov effects will alter this conclusion. 
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3.2 Including Sudakov effects 

As was shown above, it is straightforward to calculate the NLO Drell-Yan cross section ac- 
counting for the k t < k cut in the NLO contribution, and to use the scale //p = /io(^o) in 
the LO part. However, when ko <C M we have to take care of the resummation of possible 
double logarithmic terms, (a s In 2 (M / ko)) n , which account for the small probability not to emit 
additional 'bremsstrahlung' gluons with k t larger than ko- In the double Log approximation 
the corresponding Sudakov form factor reads 

T q = exp(-a s S q (fx , /i T )) (22) 

with 

^ = ^ln 2 (/i T //i ) , (23) 

7T 



where Cp = 4/3. (For an incoming gluon the corresponding expression (23) contains the colour 
factor Ca = N c = 3 instead of Cp.) The upper limit \x? of the k t integration in the T 2 factor 
is typically chosen to be \it = M. This value provides the correct single logarithmic term in 



T-factor j9]. Moreover, to be precise, expression (23) is actually of the form 



S q = ^m 2 (-/4/Vo) = — Hait/ao) + W2f • (24) 

So we include the 7r 2 -enhanced term — Cftt/A in S q . As was discussed in [TP] , this term arises 
from soft gluon emission and can be exponentiatecQ 

Certainly, the LO contribution should be multiplied by the T 2 factor given above. On the 
other hand, if one takes the complete NLO contribution, the first a s term of T-factor was 
already included in the NLO loop calculation. Therefore, to avoid double counting, we must 
compensate for this term. That is, we replace the NLO contribution C NLO by 

Cl LO -> C* LO + aASafaipd+SbfafPr)] • (25) 

Note that for the subprocess qg — > q + 'J* (that is, for a = q and b = g) we allow for the different 
probabilities of bremsstrahlung from gluon and quark lines. We choose the upper scale in the 
Sudakov T-factor to be the mass of the final parton system, that is a*t = M in the LO qq — > 7* 
case and = s for the NLO qg — > q + 7* subprocess. 

Formally working to NLO accuracy we have to keep the QCD coupling a s fixed and to use 



the double Log approximation (23). On the other hand, at relatively low scales the effect of 



the a s running may be not negligible. The effect is shown by the dashed curves in Fig. 



which are obtained by using a more precise expression (12), and the analogous expression for 
incoming gluons. This leads to a small increase in the values of /xq for k > M. Moreover, it 
is clear from Fig. [2] that in the region of interest, k in the interval (M/8, M), the possibility 
of additional QCD radiation (the Sudakov effect) does not change the value of /jlq too much. 
Recall also that the variation of the upper scale \lt corresponds to emission of a hard gluons. 
We cannot justify the exponentiation of this contribution. It should be considered as a higher 
a s order correction and should be accounted for when calculating the NNLO terms. 



4 The imaginary part is cancelled between the amplitude A and its complex conjugate A* . 
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Figure 3: The 'LO+NLO' cross sections for the Drell-Yan production of a pair of mass M — 6 

GeV and M = 12 GeV, as a function of its rapidity at the LHC energy of 8 TeV. The lower set of 
curves correspond to the imposition of a cut k t (fi + fi~) < M/4, where the optimal scale is given 
by the curves in Fig. [2] with the Sudakov effects included, namely /x = 0.42M. The upper set 
of curves correspond to the absence of a k t cutoff, in which case the optimal scale /x = 1.45M, 
see Section [2] The la error corridors obtained using the MSTW NLO parton set are shown. The 
scales in the left and right hand sides of the figure have been adjusted to allow for the naive 1/M 3 
behaviour of the cross section, so that any departure from the 'symmetry' of the plot is due to the 
PDFs. 

4 Predictions for low-mass Drell-Yan production 

The cross sections expected using MSTW08 NLO partons are plotted as the function of Y in 
Fig. [3] for M = 6 and 12 GeV for the LHC energy a/s=8 TeV. The upper curves correspond to 
the cross section obtained without the imposition of the kt cut (and correspondingly without 
the T- factor, i.e. T = 1), while the lower curves are calculated imposing a cutoff k t < M/4 on 
the transverse momentum of the produced /x + fi~ system. For the LO part the optimal scales 
/i = 1.45M (without the cut) and fx = 0.42M (with the cut) were used, as deduced in Sections 
[2] and [3] respectively. 
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Figure 4: The sensitivity of the Drell-Yan yU + /i~ cross section with respect to variations fi F = 
Ho/ 2, no in the scale of the PDFs convoluted with C^°. The cutoff k t {n + 1 1 ) < M/A is imposed, 
with optimal scale Ho = 0.42M as given by Fig. [2j 



Figs. [4] and [5] show, respectively, the sensitivity of the Drell-Yan cross section, for the 
production of a pair of mass M=6 or 12 GeV and rapidity Y, on the factorization scale 

Hf and on the choice of NLO PDFs O [TTJ, together with the error corridor obtained using the 
MSTW partons. 

In Fig. [4] we plot the cross sections with the cutoff k t < Mj 4 calculated using the scale 
Hf — Ho in the LO part and the factorization scale — A*o (solid lines) or /xp = /Jo/2 (dashed 
lines) in the remaining NLO part. Just as for the cross sections without the cutoff, shown in 
Fig. [TJ the use of the optimal scale /i in the LO part provides sufficient stability with respect to 
the variation of the factorization scale in the remaining NLO contribution. We do not consider 
here the domain of large Hf > Ho > ^0 hi the remaining NLO contribution, since then a large 
transverse momentum of the incoming parton immediately violates the cutoff condition k t < k . 



5 Conclusion 

It is well known that the Drell-Yan cross section for low-mass production depends sensi- 

tively on the choice of factorization scale, \Lf- in an earlier publication [1] we showed how this 
sensitivity can be greatly reduced by fixing the scale hf by comparing the exact contribution 
of NLO gq — > q^* subprocess with the analogous contribution generated by LO DGLAP evolu- 
tion. Since the NLO gq — > qj* subprocess is described by the same Feynman diagram as that 
in the last step of the evolution, the scale /i at which the LO evolution reproduces the exact 
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Figure 5: The sensitivity of the Drell-Yan cross section with respect to the choice of PDFs 

[HI [IT]. The cutoff k t (/-i + [i~) < M/4 is imposed, with optimal scale /xo = 0.42M as given by Fig. 
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NLO result may be treated as the factorization scale probed by the Drell-Yan process. After 
the scale /i is fixed in the LO part, the prediction is found to depend weakly on the variation 
of factorization scale in the remaining NLO part. In Section [2] we illustrated the procedure 
by calculating the value of the optimal scale analytically, and found fio = 1.45M, where M 
is the mass of the pair. Experiments have so far reached down to M ~ 6 GeV and up 

to rapidities Y ~ 4, which probe parton PDFs as low as x ~ 10 -5 , albeit at high Q 2 scales, 
Q 2 = (1.45M) 2 ~ 75 GeV 2 for M ~ 6 GeV. 

In this paper, we show that lower scales can be probed if we impose a cutoff on the transverse 
momentum of the pair, k t < k Q . In this case the optimal scale depends on the choice 

of the cutoff k . First, in Section 3.1 we performed the analytic calculation of fj, (k ), see 
Fig. [2] Then, we included the effects of the Sudakov suppression T-factor. Fig. [2] shows that 
in the region of not too low scales (corresponding to ko > M/8) the T-factor makes only a 
small increase in the value of the optimal scale. For example if we take ko = M/4, then the 
Sudakov effect only increases /j,q from 0.39M to 0.42M. Hence M ~ 6 GeV data probe PDFs 
at Q 2 ~ 6 GeV 2 . The cost of using this cutoff in shown in Fig. [3] that is, imposing the cut 
kt < M/4 GeV on the M = 6 GeV data, means we keep only about a third of the events (and 
for M = 12 GeV we keep about half of the events, but really, to reach lower scales, we should 
now choose a lower cutoff, say ko ~ M/8) 

Fig. [3] is computed using MSTW NLO partons [5], and shows the error corridors of the 
predictions. As expected the percentage uncertainty is largest at the lower Q 2 scales, which 
indicates the value of imposing the cutoff to reveal information of the low x behaviour of PDFs. 
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Fig. [5] shows the range of predictions obtained using other recent sets of PDFs [TT] . 

What is the impact of such measurements on global PDF analyses? Clearly, the measure- 
ments can probe a small x domain so far unexplored by data. If, at such small x, we were to 
believe in pure DGLAP evolution, then the majority of the low x quarks (antiquarks) comes 
from the g — )■ qq splitting. Therefore the M dependence of the Drell-Yan cross section will 
probe the gluon distribution in a similar way to how the scaling violations in deep inelastic 
scattering data, dF 2 (x, Q 2 )/dQ 2 , probe, via DGLAP evolution, the low x gluon density. Note 
that for usual fully inclusive kinematics, at a fixed initial energy, changing the value of M 
simultaneously changes the value of x. Now, by selecting events with different k t cuts, we may 
vary the scale [if = /^o while keeping the value of x fixed. 

Figs. [3]and[5]may give a misleading impression about the value of low-mass Drell-Yan data. 
Data for the production of pairs of mass M < 20 GeV at rapidities Y > 3 probed the 

unexplored domain x < 10~ 4 . In this x region the PDFs curves in Fig. |5]are based on unreliable 
extrapolations assuming pure DGLAP forms. Clearly the extrapolations must be modified to 
allow for absorptive effects and ln(l/a;) corrections, etc. Indeed, it is fair to say that there is no 
meaningful prediction in this domain. Rather we should regard the low-mass Drell-Yan data 
at the LHC as offering an important direct measurement of the PDFs in this unexplored low x 
region. 
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